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, The two-qubit canonical decomposition SU{4) = [SU {2)CX>SU (2)]A[S{7(2) (S)S{7(2)] writes any two-qubit unitary 

O , operator as a composition of a local unitary, a relative phasing of Bell states, and a second local unitary. Using 

■ Lie theory, we generalize this to an «-qubit decomposition, the concurrence canonical decomposition (CCD) 

SU{2") = KAK. The group K tixes a bilinear form related to the concurrence, and in particular any unitary in 
K preserves the tangle | (0| (— I'dj) ■ ■ ■ {—i(5),)\i)) p for n even. Thus, the CCD shows that any n-qubit unitary is a 
composition of a unitary operator preserving this n-tangle, a unitary operator in A which applies relative phases 
to a set of GHZ states, and a second unitary operator which preserves the tangle. 

As an application, we study the extent to which a large, random unitary may change concurrence. The result 
states that for a randomly chosen a £ A C Sf/(2-''), the probability that a caiTies a state of tangle to a state 
of maximum tangle approaches 1 as the even number of qubits approaches infinity. Any v = ki ak2 for such an 
a e A has the same property. Finally, although | ((|)| (— /Oj ) ■ ■ ■ (— roJ,)|i|)) p vanishes identically when the number 
of qubits is odd, we show that a more complicated CCD still exists in which K is a symplectic group. 
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I. INTRODUCTION 

+-> ' 

^ ■ Entanglement is a unique feature of quantum systems that plays a key role in quantum information proccessing. 

^ ' Much effort has gone into describing the entanglement present in the state of a quantum system composed of 

two or more measurably distinct subsystems. Because different amounts of entanglement may be shared among 
the various partitions of the tensor factors of the Hilbert state space, there is no single measure of entanglement 
that captures all non-local correlations for many-particle systems. Rather, the number of partitions of the tensor 
factors grows exponentially with the number of factors themselves. Thus, it is reasonable to guess that same 



X 



H I is true for the number of useful entanglement measures. In fact, the situation is yet more complicated. Many 

reasonable definitions create uncountably many entanglement types, which thus may not be associated to countable 
collections of partitions or monotones. 

Nevertheless, it is interesting to consider how much entanglement is created by a given unitary evolution U 
of an «-qubit state space. To achieve this in a limited context, we focus on a single multi-qubit entanglement 
measure, the «-concurrence (1). Using Lie theory, we may decompose a unitary operator acting on n qubits into a 
form such that the entangling power of the unitary with respect to this measure is manifest. 

The «-tangle and its square root, the n-concurrence, are two of several proposed multiqubit entanglement mea- 
sures. Others include polynomial invariants which involve moments of the reduced state eigenvalues (2), the 
Schmidt measure (3) which is related to the minimum number of terms in the product state expansion of a state, 
the Q measure (4) which is related to the average purity each qubit's reduced state, and GAVIN ADDS SOME- 
THING. (5) A further measure makes use of hyperdeterminants (6); this powerful technique makes computation 
difficult in more than six qubits. The concurrence C„ is originally introduced in the two-qubit case (7). It is 
generalized to a measure on two systems of arbitrarily many dimensions in {■:,) and extends to «-qubits ( ). 

We now consider the quantitative expression for the concurrence. Suppose a quantum state space of data for 
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a quantum computer. Specifically, fix « as the number of qubits, = 2". Throughout, we use \ j) not to denote 
the state of a qudit but rather as an abbreviated multi-qubit state via binary form. For example, in three qubits 
|5) = |101) = |1) (g) |0) (g) We write ^„ = spanj,{|;) ; < ; < A^- 1} for the «-qubit Hilbert state space. Then 
the concurrence is a map C„ : 9{„ [0,°°) given by C„(|\|/)) = |(\(/|(— /Oj) • • • (— /0'^,)|i(/)|. Note that the expression 
inside the complex norm is in general not real. A related entanglement measure x„ = C,^ is known as the n-tangle 
when n is even For = 1, the «-tangle of a state assumes real values in the range < T„ < 1. It is 

moreover an entanglement monotone (9), as any good measure should be. This means in particular that T„ : 9{„ 
[0,°o) vanishes on full tensor products of local states, and moreover that T„(vi (g V2C3i • • • (g v„|i|/)) = T„(|\|/)) for any 
vi (8) V2 <g) • • • (g v„ S ^"SU (2). We show in Appendix C that the «-concurrence is also an entanglement monotone. 

The n-concurrence only detects certain kinds of entanglement. Specifically, while it returns zero on all separable 
states, it may also return zero on certain non-separable states. We illustrate the monotone's behavior by exam- 
ple. First, the «-partite Greenberger-Horne-Zeilinger (GHZ) state |G//Z„) = (1/ V2)(|0i . . .0„) + | li . . . 1,,)) has 
maximal ^-concurrence while C„(|G//Z„_i) (g) |0)„) = 0. As a second example, the generalized \W) state given by 

I = ( 1 / ^/«) (|10...0) + |010...0)H h|0...01)) has zero «-concurrence despite being entangled. States with 

subglobal entanglement can also assume maximal «-concurrence; C„{\GHZ„)) ~ C„{\GHZ„i2} ® \GHZ,^i2}) — 1. 
Generally, the «-concurrence seeks out superpositions between a state and its binary bit flip. 

We extend these definitions by introducing a complex bilinear form, the concurrence form Cn '■ x ^4 — * C. 
Here, complex bilinear means the function is linear when restricted to each variable. The antisymmetric concur- 
rence form ,— ) is nonzero even in the case « is odd, although of course C2p-i = since C2p-i(|v)> l¥)) = 
-C2p-i(|v),|v))=0. 

Definition I.l The concurrence form C„ : :K„x ^„^C is given by C„(|v|/), = (\|/|(— /oj)(-;0-2) • • • (-/oJi)|(|)). 
Note the complex conjugation of the lead bra is required for complex linearity (rather than antilinearity) i n the first 
variable. The concurrence quadratic form is = (r„(|i(/), so that C„(|\(/)) = = Vt„(|\(/)). 

Note that Q^^ is a complex quadratic polynomial on the vector space itf„. 

The main technique of this paper is to build a new matrix decomposition of the Lie group of global phase normed 
quantum computations SU (N). It is optimized for the study of the concurrence and «-tangle and generalizes the 
two-qubit canonical decomposition (10; 11; 12; 13; 14; 15) 

SU{4) = [SU{2)®SU{2)]A[SU{2)®SU{2)] (1) 

Here, the commutative group A applies relative phases to a "magic basis" (7; 12; 13; 14; 16) of phase-shifted Bell 
states. This two-qubit canonical decomposition is used to the study of the entanglement capacity of two-qubit 
operations (15), to building efficient (small) circuits in two qubits (14; 17; 18), and to classify which two-qubit 
unitary operators require fewer than average multiqubit interactions (17; 18). 

The canonical decomposition is itself an example of the G ~ KAK metadecomposition theorem of Lie theory 
(19, thm8.6,§VII.8). This theorem produces a decomposition of an input semisimple Lie group G given two further 
inputs: 

• a Cartan involution (19, §X.6.3,pg.518) : g ^ g for g = Lie(G). By definition, satisfies (i) 0^ = 1 and 
(ii) Q[X,Y] = [0X, 07] for all X, F G g. We write g = p 1 for the decomposition of g into the +1 and - 1 
eigenspace of 0. 

• a commutative subalgebra a C p which is maximal commutative in p. 

Then write K = exp t,A= exp a, where for linear G C GL(«, C) the exponential coincides with the matrix power 
series on each of the Lie subalgebras i, a. The theorem asserts then that G = KAK = {kyako ; ki,k2 E K,a E A}. 

For example, the cononical decomposition of SU{4-) arises as follows. Take : 5u(4) su(4) by Q{X) = 
( - ia-i ) ( - /0-2 )X ( - ia\ ) ( - ia^ ) and 

a = spanR( /|0)(0| -/|2)(2| +/|3)(3|, /|0)(3| +/|3)(0|, /|1)(2| +/|2)(1| 1 (2) 



It is not clear how to recover the celebrated 3-tangle from this construction. By this definition we will see that z,, = for all odd n, so 
references to T„ in the present paper will suppose that is not the case. 
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We extend this particular construction to n-qubits. 

Definition 1.2 Let S = {-ia\){-ia2) ■ ■ ■ (-/O"^). Define 9 : su(2") ^ su(2") by Q{X) = S'^XS = {-l)"SXS. 
Then i denotes the +l-eigenspace of 6 while p denotes the — 1-eigenspace. Finally, in case n is even we define 

a = spanR({ /| j) (;| + /|A^ - j - 1) (A^ - ; - 1 1 - i\J + I) {] + l\ - i\N - j -2){N - j - 2\ ; Q < j < 2"-' - 2 } 
U { {N-j-l\+i\N~j-l){j\;0<j<2"-'^-l},) in case « even 

(3) 

with A = exp a. In case n odd, we drop the second set: 

a = spanR({ {j\+ i\N - j - I) {N - j - l\ - i\j + I) {j + l\ - i\N - j -2) {N - j -2\ ■ < j < 2«-' - 2 }) 
in case n odd 

(4) 

Modulo checks reserved for the body, the concurrence canonical decomposition fCCD j in «-qubits is the resulting 
matrix decomposition SU (2") ~ KAK. Note that n may be even or odd. 

In n-qubits, it is certainly not the case that K is the Lie group of local unitaries. Nonetheless, we prove momen- 
tarily by direct computation that the local unitary group S'C/ {2)®SU{2) (g) ■ ■ -^SU (2) C K, with strict containment 
in « > 3 qubits by a dimension count. Moreover, for « =2p an even number of qubits the concurrence canoni- 
cal decomposition is computable via an algorithm familiar from the two-qubit case (14), (see Appendix A). The 
following theorem provides the key to interpreting this extended canonical decomposition. 

Theorem 1.3 Let K = exp t for t the +l-eigenspace of the Cartan involution Q{X) — S^^XS. Then K is the 
symmetry group of the concurrence form C„- Specifically, for u E SU{N), 

{u eK)^[ Cn{u\<\?),uW)) = C„(|(^>,|¥>) for every |(^),|\(/) € ^n] (5) 

Moreover, for n even the concurrence form is symmetric. In the even case, it restricts to the usual dot product on 
the M-span of a collection ofn concurrence one states, and this R subspace of J{„ is preserved by K. On the other 
hand, for n odd Cn is antisymmetric, i.e. a two-form. Thus, 

• K = Sp{N/2), if n is an odd number of qubits 

• K ^ SO{N), ifn is an even number of qubits 

Remark 1.4 Bremner et al (20, thmS) observe symplectic Lie algebras independently in a context related to the 
above. We explore this in more detail in a sequel manuscript. O 

This interpretation allows for an extension of prior work on the entangling capacity of two-qubit unitaries ( 1 5). 
Here is the precise result: 

Definition 1.5 The concurrence capacity of a given «-qubit unitary operator v e SU (N) is defined by k(v) — 
max{C„(v|v)) ; C„(|v)) =0,(v|v) = 1}. 

Corollary 1.6 ( of 1.3) Let u = kiak2 be the n-qubit canonical decomposition of u G SU{N). Then k(m) ~ K{a). 

Given the CCD, the function K is properly viewed as a function on the A factor rather than on the entire group 
of phase-normalized unitaries SU [N). Finally, a careful analysis of K(a) for randomly chosen a in A produces the 
following, perhaps surprising result. 

Theorem 1.7 Suppose the number of qubits is even, i.e. n = 2p. Then for large p almost all a € A have maximal 
concurrence capacity. Specifically, suppose we choose a A at random per the probability density function given 
by the unit normalized Haar measure da. Then 

limp^aoProbability[K{a) = I] = limp^oo da{{a Cz A ; K{a) = I}) = 1 (6) 

We rephrase this result colloquially. Suppose we think of those states in even qubits with T„(|\(/)) = 1 as 
GHZ -like. Then as the even number of qubits grows large, almost every unitary evolution will be able to produce 
such a maximally concurrent GHZ -like state from some input state of concurrence. 
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Notation and Contents 

We provide some samples of our notation for the reader's convenience. Throughout, n is a number of qubits 
and = 2". For v = L^^io^^.'^l'^) (-'I' have the adjoint = Y!^jkl=Q^k.j\k){j\. We also require the transpose 
operation, most easily visualized in matrix form as [vjj^Y ~ i'^kj)- Equivalently, = Y!jklQ'^j,k\J)i^\- Thus 

= v^. Recall also the convention of collapsing the binary for an integer inside the ket of a computational basis 
state. We use lower rather than upper case letters for most operators to avoid confusing them with Lie groups 
denoted by capital letters. The older term scholium is used to refer to a corollary of the proof of a theorem or 
proposition rather than its formal statement. Besides these conventions, we follow the notations of either (2 I ) or 
(19). 

The paper is structured as follows. In §11, we verify that the conventions are Definition 1.2 are appropriate for 
invoking the G = KAK theorem. Having verified that the matrix decomposition exists, §11 further describes entan- 
glers and finagle rs, loosely similarity matrices which rotate the CCD onto more standard KAK decompositions of 
SU (N). In §111, we discuss the concurrence capacity and prove the properties of this capacity asserted above. The 
three appendices consecutively (i) provide an algorithm for computing the CCD given a matrix v e SU {N), exclu- 
sively in the case n is even, (ii) argue that any two normalized states |(|)), \\\t) with identical concurrence must have 
k\(^) = for some k in the symmetry group K, and (iii) prove that the concurrence C„(— ) is an entanglement 
monotone. 



II. ENTANGLERS, FINAGLERS, AND CARTAN INVOLUTIONS OF su(A') 

This section has two goals. First, we show our KAK decomposition is well-defined, by noting that 9 is a 
Cartan involution, checking by direct computation that a is abelian, and arguing that a is maximal commutative. 
Second, we prove Theorem 1.3. There are generally two approaches to the theorem. We could recall standard 
Cartan involutions and KAK decompositions from the literature. We will shortly construct similarity matrices 
Eq and Fq which rotate the standard G = KAK decompositions of SU{N) onto the CCD, and we could simply 
appeal to these matrices and the standard structures. Alternately, (many) intrinsic computations would suffice to 
check the required properties for G = KAK. The present approach is a compromise. The argument that the CCD 
SU{N) = KAK is well-defined is intrinsic, except for a single appeal to classification. On the other hand, the 
classification of the K groups uses similarity matrices. As such, it is ultimately a change of basis in the n-qubit 
state space 



Properties of the ccd SU{N) = KAK 

The following proposition is not used in the sequel. However, we include a direct proof due to its importance 
in guiding the choice of 0. It simplifies an older argument and arose from correspondence with P.Zanardi. 

Proposition II.l Let K be as in Definition 1.2. Then there is an inclusion SU (2) ®SU{2) ® ■ ■ ■®SU{2) C K. 

Proof: Recall i(f = /|0)(1| +/|1)(0|, /o-^' = |0)(1| - |1)(0|, and i(5~ = /|0)(0| - forms a basis of su(2'). 
For the statement of the proposition, it suffices to check Lie[(E)i5't/(2)] = spanj/CTy , /a^, /o^ ; 1 < j < «} C t. We 
further recall the last item of Lemma II. 2, as well as the fact that the complex conjugates of the Pauli matrices are 
i(5^ — ~i(f, ia^ — i& , and ia' = —i&. Then we wish to show that fixes every cr^, o^, and o^. For this, 

(-1)" 5(g) 5 = (-l)«5Hop5 = i-irsHia-]) = (ia-)) 
i-iy Siia^S = (-l)«5(/o^)5 = (-l)"52(/op = (/op 
(-l)"5(/o))5 = i-irS{-i&j)S = i-iysHio]) = (ia)) 

Hence each such infinitesimal unitary is in the +1 eigenspace of 0. This concludes the proof. □ 
We next note that is a Cartan involution. Indeed, direct computation shows that 0^ = 1. Moreover, 



[QX,QY] = {S-^XS)iS-^YS) - {S-^YS)iS-^XS) = S-\XY -YX)S ^Q[X,Y]. (7) 
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Thus we need the following to complete the argument that SU (N) = KAK of Definition 1.2 is well-defined: (i) 
a C p, (ii) a is commutative, and (iii) no larger subalgebra of p containing a is commutative. 

Lemma II.2 Ler5= (-/aj)(-/a2) • • • (-/oJi) be as in Definition 1.2. Then (i) S\j) = (-1)*-'|A'-; - 1), (ii) {j\S^ 
(— 1 )"^*-' (A^ — J — 1 1, and ( iii) So'j = — 0'y5, S&j = 0^5, and S&j = —o'jS. Note that ( ii) refers to a composition of 
linear maps. 

Sketch: For (i), compute. For (ii), consider {j\S\k) for \k) varying over all computational basis states. Then apply 
(i) for (ii). For (iii), nonlike Pauli matrices anticommute, while S itself is a tensor □ 

Lemma II.3 Let a be as in Definition 1.2. Then a C p. 

Proof: There are two coordinate computations to complete in this case. For the first, momentarily extend the 
definition of 9 to 9 acting on u(A^) by the same formula. Then 

9[ i\j) 01 + /|A^ - j - 1) (A^ - - 1 1 ] = {-rrS{ -i\j) {j\ -i\N-j-l){N-j-l\)S 

= {-ir'i[{-ir\N-j-i){N-j-i\+i~i)"\j){j\] 

= -i\j){j\-i\N-j-l){N-j-l\ 

Thus i\j) {j\ + i\N — j — ^) {N — j — l\ is in the — 1 -eigenspace of 9, so that the elements of the first set of the 
definition of a are contained in p. For the second basis set in case n even, 

9[ /|;) (A^ - - 1 1 + /|A^- - 1) (j| ] = 5[ (-/)|;) (A^ - y - 1 1 + (-/) |A^ - y - 1) 01 ]S 

= [ i-l)*-'+^"-^"-^J^\-i)\N - j - l){j\ 

+ {-l)"-*j+^"-*j\-i)\j){N - j - l\] 

= H)|A^-y-i)01 + H)0)(A^-y-i| 

Thus i\j) (A' - J - 1 1 + i\N- j - 1) 01 G P> in case n even. □ 
Proposition II.4 Recall a from Definition 1.2. Then a is commutative. 

Proof: Throughout, < j,k < N/2 — l,N = 2". The following three computations of Lie brackets suffice. 

[ /|y>OI+'lA^-i-l>(A^-y-l| , i\k){k\ + i\N-k-l){N-k-l\ ] = 
-\j){j\k){k\-\N-j-\){N-j-l\N-k-l){N-k-l\ + 
\k){k\j){j\ + \N-k-l){N-k-l\N-j-l){N-j-l\ 



[ Hr^'\J){N-j-l\+i"-'\N-j-l){j\ , {~i)"+'\k){N-k-l\ + i"-'\N-k-l){k\ ] = 
-\j){N-j-l\N-k-l){k\ -\N-j-l){j\k)(k\ 

+ \k){N-k-l\N-j-l){j\ +\N-k-l){k\j)(N-j-l\ 



[ i\j){j\+i\N-j-l){N-j-l\ , {-ir+'\k){N-k-l\ + i"-'\N-k-l){k\ ] = 
{-i)"\j){j\k){N-k-l\ + i"\N-j-l){N-j-l\N-k-l){k\ 
-{-iY\k){N-k-l\N-j-l){N-j-l\ - i"\N-k-l)(k\j){j\ 

Each of the final expressions is zero in case j k and also zero in case j = k. Thus, a is commutative. □ 

The arguments above almost complete the proof that the CCD SU (N) = KAK is well-defined. In the abstract, 
one also needs a fairly large coordinate computation which verifies a is maximal commutative. This would verify 
that for any X G p with [X,H] =0 for all H E a, one must in fact have X E a. 

Rather than complete that task, we instead appeal to the Cartan classification ( , pg.5 18, table V). Ostensibly a 
classification of globally symmetric spaces, this classification also describes all possible Cartan involutions of any 
real semisimple group G up to Lie isomorphism. For G = SU (N), there are three overall possibilities grouped as 
type AI, All, and AIII. For each, the rank refers to the dimension of any maximal commutative subalgebra a of 
p. This dimension may not vary by subalgebra, since any two such ai, a2 must have kaik^^ = 02 for some k E K. 
We now excerpt from the table the possibilities for G = SU (N): 
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type 


domain g of 9 : g ^ g 


isomorphism representative of K 


rank 


AI 


su(A^) 


SO{N) 


A'- 1 


All 


su(A^) 


Sp{N/2) 


N/2-1 


Ain 


su(A^) 


S[U{p)®U{q)],p + q^N 


mm{p,q) 



Suppose then for the moment that the number of qubits « is even. No type AIII Cartan involution admits an a of 
dimension N —I. Indeed, if p + q — N, then min(;5, ^) <N/2 <N — I. The same is true of type All involutions, 
i.e. N — I > N/2 — 1. Hence we see that A must be maximal, and for « even the Cartan involution must have 
type AI. 

What remains is to prove that a is maximal in p in case n odd. This follows by a dimension count ;/ the Cartan 
involution is type All. We thus postpone noting this point until after the proof of Theorem 1.3. See Remark II. 19. 

As an aside, type AIII involutions do not appear in this work but have been used in quantum circuit design. 
Indeed, the CS -decomposition (22; 23) is an example of a KAK decomposition arising from a type AIII involution. 
Elements within the appropriate K group may be interpreted as products of computations on the last n — 1 lines 
with computations on these lines controlled on the first qubit. 



Entanglers 



In the two-qubit case, the following computation E has the following property: 

/ 1 / \ 





Vl 



1 i 
-1 i 
-i 0/ 



satisfies E''[SU{2)(g}SU{2)]E = 50(4) 



(8) 



Using more Lie theory terminology, recall the adjoint representation of G on g given by Ad(§)[X] = gXg^^. 
Then we may restate {Ad{E^)}[SU (2) (E)SU{2)] = 5(9(4). This provides a physical interpretation for the low 
dimensional isomorphism su(2) ©su(2) = so(4). We would like entanglers for the concurrence canonical decom- 
position. 

Definition II.5 Let 9ai : su(2") su(2") denote the usual type AI Cartan involution 9ai(^) = X associated to 
SO{N) C SU (N). We say E G SU (2") is an entangler iff the following diagram commutes: 



5u{N) su{N) 
Ad{E) i i Ad(£) 

su(A^) su(A^) 



(9) 



In particular as both groups are connected, we must have Ad{E)[SO{N)] = K. 

We next prove the surprising fact that there are no entanglers when n is odd. For this, we need to recall the 
central subgroup Z [SU (N)] — {v E SU (N) ; vuv^ = u for all u G SU (N)}. The center is in fact the set of all phase 
computations corresponding to the A^*-^ roots of unity: 



Z[5C/(A^)] = {^1; = 1} (10,pg.310,516) 
With this fact recalled from the literature, we have the following lemma. 



(10) 



Lemma II.6 Suppose that for v G SU{N), [Ad{v)]{X) — vXv^ —X for every X G su(A^). Then v = some 
^ G C with = L (Hence ^ = e^'^'V^', <k<N-l.) 

Proof: Recall that exp : su{N) SU (N) is onto. Thus each u G SU (N) may be written as expX for some X. Thus, 
consider the one-parameter-subgroup ( 1 pg. 104) y : R — > SU (N) given by f i— > v[exp(fX)] v^. This has derivative 
^|f=o = vXv^ ~ X, and by uniqueness of one-parameter-subgroups (19, pg.l03,Corl.5) vexp{tX)v^ — exp{tX) 
for all t. Taking t — 1, we see vuv^ = m for a generic u G SU{N). □ 
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Proposition II.7 If the number ofqubits n is odd, then there does not exist an entangler E G U{N). 

Proof: Assume by way of contradiction that there does exist an entangler E for n odd. Then for all X e su(A^), 
we have the following equation. 

{EE^)X{EE^y = EQm[E^XE]e'^ = Q{X) = SXS'^ (11) 

Since we may vary Y = X over su(A^) as well, this implies that S^^EE^ satisfies the hypothesis of Lemma 11.6. 
Thus S^^EE^ = ^1 for = 1 or EE^ = (^l)'^- Contradiction, for EE^ is always a complex symmetric matrix 
while (^1)5 is not a complex symmetric matrix when n is odd. □ 

Scholium II.8 For an even number of qubits n, the matrix E G U{N) is an entangler ijf EE^ = (^1)'^> where 

There are many possible entanglers. Indeed, even in two-qubits other choices have been used (10; 15). One 
possibility given n even is to take the n/2 fold tensor product E ®E ® ■ ■ ■ ®E . However, we prefer the following 
choice as a standard instead, since it highlights the mapping of the computational basis to Greenberger-Horne- 
Zeilinger states. 

Definition II.9 Suppose n is even, and write S = (— /o-[)(— /Oj) • • • (—/On) = ifjlo^^ {N - j - 1| + |A^— j- 
1) with Ej = (—1)*', where #j is the number of I's in the binary expression for j. The standard entangler Eq 
in «-qubits is then given by 

1 N/2-1 

^0 = -- L l;>(2;| + /|y)(2; + i|+s,(|A'-;-i>(2y1-'lA'-y-i>(2; + i|) (12) 
V2 ^=0 



Proposition 11.10 The standard entangler Eq is an entangler. 



Proof: First, we omit due to reasons of space a set of row operations which verifies that det(£'o) = 1- Then we 

■T 




may write out an expression for E^ by reversing the indices in each bra-ket pair: 



I iV/2-l 

eI = -= \'^k){k\ + i\2k+\){k\+Zki\2k){N-k-\\-i\2k+\){N-k-\\) (13) 
v2 <.=o 

Then Scholium II. 8 shows that the following computation suffices to prove that E is an entangler. 

E,El = \Y.%V\j){j\+^j\j){N-j-l\ + i'\j){j\+^j\j){N-j-\\ 
+ Zj\N-j-\){j\+E)\N-j-\){N-j-\\ 

+ tj{\N-j-l)U\+ i\]\N -j-\){N-j-l\) (14) 
= Z%''^j{\j){N - j-l\ + \N- 

This concludes the coordinate computation. □ 
In the next section, we will also make use of the following lemma. The computation is similar. 

Lemma 11.11 E^Eq is diagonal and real. In fact, E^Eq = |0)(0| - + |2)(2| - |3)(3| + --- 

Proof: Computing the reversed product: 

E^Eo = ^I^io"'|2;)(2;| + /|2; + 1)(2;| + /|2j)(2; + l| - |2; + 1) (2j + 1 1 

+ Epj) {2j\ - iepj) {2j + 1 1 - iej\2j + 1) {2j\ + ihj\2j + l){2j + 1 1 (15) 
- lS"'|2;)(2;|-|2; + l)(2y + l| 
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This concludes the proof. 

Example 11.12 Although this example is large, we explicitly describe the standard four-qubit entangler. 



□ 



£o = (l/V2) 
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maps to a relative phase of a GHZ state. 
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Finaglers 

There do not exist entanglers when the number of qubits n is odd, because t = 5p{N/2) rather than 5o{N). Cf. 
the as yet unproven Theorem 1.3. Yet the fairly abstract embedding K of Sp{N/2) into SU{N) might be made 
more standard. This is indeed possible, and we call the any matrix which rotates K to the standard Sp{N/2) a 
finagler 

Definition 11.13 Let 9aii : su(A^) su{N) be the standard Cartan involution (19, pg.445) fixing sp{N/2), 
i.e. eAii(^) = (g) 1n/2)X'^ {-io^ ® In/i) = i-i<^-' ® ^N/2y^X{~i&' In/i)- Then a finagler F is any 

F G SU{2") which causes the following diagram to commute; 

su{N) ^ su(A^) 
Ad{F) i i Ad{F) (17) 

su(A^) su(A^) 

If F G SU (N), then we say F finagles iff F is a finagler. 

Proposition 11.14 F is a finagler ijfF{-i&- ® In/iVf'^ = ■ ■ ■ {-io-,) = (^1)5, t,^ = 1. 

Proof: For convenience, label E = -io^ ® l^/j. (F finagles) [FT.-'^ (fV(F)I.F'^ = S'^XS VZ e su(A^)] 
[FTJf^ = {t,l)S,t,^ = I]. Note that the second equivalence uses Lemma 11.6. □ 
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Example 11.15 In three qubits, we see the following computation is a finagler by direct computation. 



(1/V2) 



/ 1 


1 








o\ 


10 





-1 








10 








-1 





1 











1 


1 











-1 


10 








1 





10 





1 








V 1 


-1 








0^ 



(18) 



Unlike entanglers, it is possible for F = F. The finagler maps computational basis states to GHZ states. O 

Definition 11.16 Fix n an odd number of qubits. Let S ~ (-/crj)(— /O2) • • • (-'0>i) = L^io ' ^jiW^ J ^ ^)U\ ~ 
\j) {N — j — 1 1) with ij = ±1. The standard finagler Fq is defined to be the following linear operator: 



N/2-1 

Fo= E |;)01 + |A^-j-l)01 +l,(|j)(A^/2 + jHA^-j-l)(A^/2 + j|) 



(19) 



Note that the standard finagler is real. 

Proposition 11.17 The standard finagler Fq finagles. 

Proof: We again omit the column operations verifying det(Fo) = 1, as this would take several pages. Thus, let 
E = — ;a-*' (8) In/2 expanded as E = ^Jio ' \J) (-^/^ + J\ ^ W/2 + j) {j\. We have the following equation: 

I N/2-1 

FqL=— E |y)(A^/2 + ;| + |A^-;-l)(A^/2 + ;i -l,(|j)(j|-|A^-y-l>01) (20) 
v2 j=Q 



Moreover, F^ = ijZ^ ' |;) {j\ + \j) {N-j-\\+ ij{\N/2 + j) {j\ - \N /2 + ;) (A^- y - 1 1). Thus we see that 



\i:%''ij{\j)U\-\j){N-j-l\ + \N-j-l)U\-\N-j-l){N-j-l\ 

- 01 + li) (A^ - i - 1 1 - |A^- 7 - 1) 01 - |A^- J - 1) (A'- 7 - 1 1) 
Z%''ij{\N-j-l){j\-\j){N-j-\\) 



(21) 



This concludes the proof. □ 
We also briefly review how Sp{N/2) embeds into SU{N). By one standard definition of the group (19, pg.446). 



sp(A^/2) 



Xi X2 

Xi -xj 



; Xj ~Xj,X2,3 symmetric 



(22) 



Another standard definition (24, pp. 34-36) uses a symmetry in matrices of quaternions. Note that the matrices of 
Equation 22 are not elements of su(A^). Rather, the +1 eigenspace of 9aii(^) = E^'XE is: 



5piN/2) 



V W 

-W^ V 



; V e u{N/2),W = is complex symmetric 



(23) 



(For example, Sp{4) C SU{^) this is 36 dimensional. For W includes two real symmetric matrices with 10 di- 
mensions each, while u(4) is 16 dimensional.) One may verify this is also a copy of sp{N/2), so that Hs a copy 
of sp{N/2) as well. Also, note that for k £ K, in particular k E ®"SU (2), we expect FkF' to be in the copy of 
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Sp{N/2) above rather than to be a real matrix in an orthogonal subgroup of SU (N). Finally, note that the expo- 
nentiating the Lie algebra above is not the best way to write out a closed form for elements of the global group 
Sp{N/2). Rather, we have a block form: 

A b\ is symmetric, B'^D is symmetric, I 



Sp{N/2) = {V e SU{N) ■V'I.V = Y.} = 



K is the symmetry group of the concurrence form 

We are now in a position to provide the physical interpretation of K. Namely, K is the symmetry group of the 
concurrence bilinear form, as stated in Theorem 1.3. 

Proof of Theorem 1.3: We first prove that v e iff (:„(v|(j)),v|\j/)) = Cn{\<^),\^)) for all |(j)), |\|/) e LetX = 
log V. Since X e su(A^), X is anti-Hermitian, i.e. X ~ —X'' — ~X^. Finally, recall S = (— ;aj)(— /a,) • • • {~ial). 
Thus in mathematical notation, we have for w,x e y{„ the concurrence form given by C„{w,x) = w^Sv. Hence 

{X = S- ^XS) ^ {SX =XS) ^ {SX = -X'^S) ^{X'^S + SX = 0)^ {v^Sv = S) (25) 

Now the first item is equivalent to v E K while the last is equivalent to c{vw,vx) = {w^v^)S{vx) = w^{v^Sv)x = 
w^Sx = C{w,x) for all w,x E 9<„. 

We next prove that for n odd, K = Sp{N/2). To do so, it suffices to show n odd implies C„ is a nondegenerate 
two-form on y{„. We first show Cn{x, w) = —Cn(w,x) for any w,x E y^,,. Noting that the transpose of a 1 x 1 matrix 
is again the same matrix, we realize that C„ is a two-form as follows: 

C„iw,x) = w'^Sx = [w'^Sxf = x'^S'^w = -x'^Sw = ~C„{x,w) (26) 

Moreover, consider the tensor expression for S. We see that no eigenvalues of S are zero, and hence the form is 
nondegenerate. Thus, we must have K = Sp{N/2). 

Suppose now that n is even. We finally prove K ^ SO{N). It suffices to construct a real vector space Vr C 
so that the following properties hold: 

• The restriction of C„ to Vr x Vr is the usual dot product in the coordinates of a given basis. 

Consider then Vr = spanRjisoly) ; < y < — 1}, for £0 the standard entangler of Definition 11.9. Since £0 is 
an entangler, certainly • Vr C Vr since K acts on this real vector space by (real) orthogonal maps. Moreover, 
consider the concurrence on Vr. For w,x in the M span of the computational basis, we have Eqw, Eqx generic 
vectors in Vr. Then 

Cn{EQW,EQx) ~ {EqwY S{Eqx) ~ w^EqSEqx = w^EqEqEqEox ~ w^lx — wx (27) 

with the fourth equality by Lemma II. 11 . Hence in an even number of qubits, K fixes a real inner product on a real 
vector subspace of Thus K = SO{N). □ 

Scholium 11.18 For n even, for Eq the standard entangler of Definition II. 9, for any |(|)), \\\f) E ^n, we have 
C„(£o|c^>,£o|¥»=M¥>- 

Remark 11.19 Note that independent of any discussion of the algebra a in n an odd number of qubits, we have 
shown that the Cartan involution 9 has type AH. Hence any commutative N/2 — I dimensional subalgebra of p 
must be maximal, and the concurrence canonical decomposition SU (N) ~ KAK is well-defined for « odd. O 

Remark 11.20 Similar to Scholium 11.18, note that the standard (real) finagler Fq of Definition 11.16 translates 
between the concurrence and the more standard two-form (w,x) i— > w^[{~io^') ® 1a?/2]jc- Indeed, Fq"' = Fq since 
Fq is orthogonal. Moreover, let w, x be in the real span of the computational basis states {\j) ; < y < A^ — 1}. 
Then we may view {FqIJ) ; Q < j < N — 1} as a finagled basis, and the pullback of the concurrence from the 
finagled to the computational basis is the model two-form. Indeed, labelling E = (— /o* ) ® Iat/i, FqLFq — S and 
Fq real imply E = FqSFq. Hence Cn{Fow,Fox) = {FowYS{Fqx) = w^Ejc. O 
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Cartan Involution in Coordinates 

We finally present the Cartan involution in coordinates and provide some sample calculations. Let X G su(A^), 
say with X = T!jkloXj^k\k){j\. We now compute explicitly Q{X) so as to arrive at coefficient expressions for p, t. 

m = i-i)"Li^ioxj,kS\k){j\s ^i-ir+"-''^'-'LUo^j,k\N-k-i){N-j-i\ 

Consequently, we have the following characterizations: 

• X = l^l-,^loXj,k\k){j\ e p iff [{X e su{N)) and (x^-i-^.^-i-,- = (-1)#-Wii^j)] 

• X = LlkloXj.k\k){j\ e i iff [{X e su(A?)) and (x^-i-^.^v-i-,- = 
This moreover produces the following description of i. 

« = span^ {\k){j\-\j){k\ + i-lfj+^'\N-k-l){N-j-l\-{-l)^J+^'\N-J-l){N-k-l\} 

U {i\k) {j\ + i\j) {k\ + (-l)#^+#>+i/|A?~ y -l){N~k-l\ + {-iyj+^''+U\N-k -l){N-j-l\} 
U{i\J)(j\-i\N-j-l){N-j-l\} 

(28) 

Remark 11.21 We warn the reader that the above expression does not allow one to count dimensions. Several 
repetitions occur from set to set, and moreover the expressions may vanish in case j = N — k — I. O 



Example in the two-qubit case 

Recall the subalgebra [1 ® su(2)] © [su(2) ® 1] of infinitesimal transformations by SU(2) (g) SU{2) C SU{4-). We 
show how the above Equation 28 recovers this subalgebra in the case of n = 2 qubits. 

We begin by plugging k = Q, j = 1. Expanding into binary (or writing out the matrix) makes clear this is a 
tensor, and moreover a tensor by an identity matrix. Recall again that both are required to be in the Lie algebra of 

SU{2)(x)SU{2). 

|0)(1|-|1)(0|-|3)(2| + |2)(3| 

|00)(01|-|01)(00|-|11)(10| + |10)(11| = (29) 
(|0)(0| + |1)(1|)®(|0)(1|-|1)(0|) 

One may similarly analyze the following matrices: 

/|0)(1|+/|1)(0|+,13)(2|+/|2)(3| 

|0)(2|-|2)(0|-|3)(1| + |1)(3| (30) 
.■|0)(2|+/|2)(0|+/|3)(1|+/|1)(3| 

Note that for the next four expressions, substitution returns a matrix: 



|0)(3|-|3)(0| + |3)(0|-|0)(3| 
/|0)(3|+/|3)(0|-/|3)(0|-/|0)(3| 
|1)(2|-|2)(1| + |2)(1|-|1)(2| 
/|1)(2|+/|2)(1|-/|2)(1|-/|1)(2| 



(31) 



Further substitution yields the following: 



|1)(3|-|3)(1| + |2)(0|-|0)(2| 
/|1)(3|+/|3)(1|+/|2)(0|+/|0)(2| 
|2)(3|-|3)(2|-|1)(0| + |0)(1| 
/|2)(3|+/|3)(2|+/|1)(0|+/|0)(1| 



(32) 
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Finally, we consider the diagonal matrices in t: 

ilO)(0|-/|3)(3| 
/|1)(1|-/|2)(2| 

Note that the R span of these two matrices coincides with R(/oJ) ® M(;o?). 

The Cartan involution formalism thus works, although in a cumbersome way. We next explore the answer it 
returns in the three-qubit case. 



Example in the three-qubit case, K = 5p(4) 

We now describe explicitly the output of Equation 28 in three qubits. The corresponding real Lie algebra is 
thirty-six dimensional, which implies by the Cartan classification that K is an abstract copy of Sp{4). A copy of 
5(9(8) would rather be twenty-eight dimensional. 

The simplest way to organize the three qubit computation is to appeal to separation. We say a term \k){j\ has 
separation \k — j\ and extend linearly. In Equation 28, each matrix described has a well-defined separation. 

Separation /|0)(0| -;|7)(7| 
Total 4 /|1)(1|-/|6)(6| 

/|2)(2|-/|5)(5| 

/|3)(3|-/|4)(4| 



Separation 1 


|o)(i|- 


|i)(o|- 


|7)(6| + 


|6)(7| 


Total 8 


/|0)(1|4 


/|1)(0| 


+ '17)(6| 


+ '16)(7| 




|1)(2|- 


|2)(1| + 


|6)(5|- 


|5)(6| 




'11)(2|4 


'12) (1| 


-'16)(5| 


-/|5)(6| 




|2)(3|- 


|3)(2|- 


|5)(4| + 


|4)(5| 




i|2)(3|4 


/|3)(2| 


+ '15)(4| 


+ '14)(5| 




|3>(4|- 


|4)(3| 








'13>(4|4 


'14)(3| 






Separation 2 


|0)(2|- 


|2)(0|- 


|7)(5| + 


|5)(7| 


Total 6 


/|0)(2|4 


/|2)(0| 


+ '17)(5| 


+ '15)(7| 




|1>(3|- 


|3)(1|- 


|6)(4| + 


|4)(6| 




'11>(3|4 


'13)(1| 


+ '16)(4| 


+ '14)(6| 




|2)(4|- 


|4)(2| + 


|5)(3|- 


|3)(5| 




'12)(4|4 


'14) (2| 


~'15)(3| 


-'13)(5| 


Separation 3 


|0>(3|- 


|3)(0| + 


|7)(4|- 


|4)(7| 


Total 6 


/|0)(3|4 


'13) (0| 


-'17)(4| 


-'14)(7| 




|1)(4|- 


|4)(1| + 


|6)(3|- 


|3)(6| 






-'14)(1| 


-'16)(3| 


-/|3)(6| 




|2)(5|- 


|5)(2| 








'12) (5 K 


-'15)(2| 







(35) 



(36) 



(37) 
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Separation 4 |0) (4| - |4) (0| - |7) (3| + |3) (7| 
Total 4 /|0)(4| + /|4)(0| + /|7)(3| + /|3)(7| 

|1)(5|-|5)(1|-|6)(2| + |2)(6| 

/|1)(5| + /|5)(1| + /|6)(2| + /|2)(6| 

Separation 5 |0)(5| - |5)(0| + |7)(2| - |2)(7| 
TotaU /|0)(5|+/|5)(0|-i|7)(2|-/|2)(7| 

|1)(6|-|6)(1| 

/|1)(6| + /|6)(1| 

Separation 6 |0)(6| - |6)(0| + |7)(1| - |1)(7| 
Total 2 /|0)(6| + /|6)(0|-/|7)(1|-/|1)(7| 

Separation 7 |0) (7| - |7) (0| 
Total 2 ;|0)(7| + /|7)(0| 

Thus we see a total of 4 + 8 + 6 + 6 + 4 + 4 + 2 + 2 = 36 real dimensions in t. Now by the Cartan classification (19, 
pg.518), the Cartan involution 9 must be either type AI fixing an abstract copy of 5(9(28), type AIII fixing some 
S[U{p) U{q)] for p + q = 2", OT else type All fixing an abstract copy of Sp{4). Since only Sp{4) is thirty-six 
dimensional, we see ? = sp(4) and K = Sp{4-). 

III. APPLICATIONS TO CONCURRENCE CAPACITY 

This section focuses on an application of the concurrence canonical decomposition SU (N) = KAK of Definition 
1.2 when the number of qubits n is even. Namely, we study how a given computation v G SU (N) may change the 
concurrence of the quantum data state. Since we have the concurrence C„(|x|/)) = | (\(/|(— /Oj)(— /O'j) • • • (— /o>,)|\|/) | 
with the «-tangle T,, = C^, for « even, there are immediate applications to the «-tangle as well. 

Let V e SU {N). Recall from Definition 1.5 that the concurrence capacity is defined as 

K(v) = max{C„(v|v)) ; C„(|v)) - 0, (viv) - 1} (42) 

Since we vary over all C„(|\(/)) = 0, we see that for k E K we have K{vk) ~ K(y) by symmetry. Immediately 
K{kv) = k(v). Thus, for v = kiak2 the C . C . decomposition of any v e SU (N), we have K(y) = K{kiak2) = K{a). 

We next describe the concurrence capacity of any a E A. The formalism makes strong use of entanglers to 
translate between d and {w,x) w^x. 

Definition III.l The concurrence spectrum 'kciv) of v G SU{N) is the spectrum of EqvEo{eIvEoY , for £0 the 
standard entangler of Definition II. 9. Note that the spectrum is the set of eigenvalues since y{„ is finite dimensional. 
The convex hull CH[Xj.(v)] of Xc{v) is the set of all line segments joining all points of Xc{v), i.e. 

cH[X,(v)] = <^ fj^j ; < < 1 , L - 1 (43) 

I zjeX,iv) j=o ) 

These definitions allow us then to prove the following general results regarding concurrence capacity. The 
techniques closely follow those in prior work ( ). 

Lemma III.2 Let v eSU (N), with CCD v = k\ak2for a = EodE^for d diagonal in SU{N). 
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. K(y) = K(«) = mflxj I I^^Jo' 02^/21 . ^ ^^-^1 aj\J),{^v\^v) = hWW) = oj 

• (k(v) = K(fl) = 1) ^ (0 e CH[?t,(v)] - CH[X,(fl)]). 
Proof: For the first item, recall v = k\ak2- Thus the following expression results from expanding Definition III.l. 

eIvEq{eIvEoY = [(4fci£o)(£(ifl£o)(£^^2£o)][(£(ife£o)^(£({fl£o)^(£(l^i^o)^] (44) 

Label the elements of SO{N) by oi = (/Sq^i/sq), 02 = (£'qA:2-£'o), and put d = EqoEq diagonal. Then the above 
reduces to o\do202d^ o\ ~ oid^o^^, with spectrum identical to d. 

For the next item, compare the two-qubit case (15, Eq.(41)) and recall Scholium 11.18. Suppose C„(|(p)) = 
0. Then per SchoHum 11.18, for |\|/) = £019) we have we see = C«(£'o£'()|(p),£'o£'dl'p)) = (vl¥)- Now for 
K(a), take |\|/) — Eq\(p). We then maximize over expressions Cn(a|9),a|(p)) = Cn{EoEQaEo\\\t) ,EoEQaEo\\\t)) = 

C„iEod\\\f),Eod\\v)) ^iMd^W)- 

The final item makes use of the Schwarz inequality. For should the concurrence capacity be maximal, there is 
by compactness of the set of normalized kets some normalized with C„(|(|))) = 1. For |\|/) = Y!j=o 

N-l N-l N-l 

i = iL«]4i^Li«'4i = EK-i' = i (45) 

j=0 j=0 j=0 

The Schwarz equality further requires some z G C, zz = 1, so that a^jdj = \aj\-^z,yj. Now since (\|/|\(/) = 0, 

N-l N-l 

0=L«'=LKI'4^ (46) 

)=o j=o 

Multiplying through by z and taking the complex conjugate, we see £ CH[A,t (v)]. □ 

As already noted in the introduction, the concurrence capacity K is properly thought of as a function of A rather 
than a function of SU (N). This is advantageous from a computational standpoint, because in order to calculate 
k(v) one need minimize over a function involving A^— 1 real parameters in A versus — l parameters describing 
a general v £ SU{N). We next consider typical values for a large number of qubits. To do so, we need to be able 
to randomly choose an element of A. 

Definition III.3 Consider the following coordinate map on the commutative group A: 

N-2 

[0,2<-i ^Aby (fo,f2,--- ,fw-2)^exp£o( L itj\j){j\- itj\j + 1\)eI (47) 

j=o 

The Haar measure on da is the group multiplication invariant measure da — (27t)^^+'iifo dt2 ■ ■ ■ dtN-2- This is 
the pushforward of the independent product of uniform measures dtj/{2K) on each [0, 2k] . 

Recall that for p = 2«, Theorem 1.7 asserts that according to da, almost all a S A have K(fl) ~ 1 for p large. 
Specifically, we assert 

limp^^da{{a G A ; K(a) = 1}) = 1 (48) 
We prove this assertion shortly, but we first need a lemma. 

Lemma III.4 Label as uniform distribution on the circle a distribution whose pullback to [0,27t] under 1 1— s- e^iw 
is uniform, and similarly say two random variables Zi, Z2 on {zz =1} are independent iff their pullbacks to 
[0,27t] X [0,27t] are. Then suppose Z\ is any random variable on the circle, and let Z2 be independent to Z\ and 
uniform. Then Z1Z2 is uniform. 
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Proof: Consider the random variable T = — /logZi — /logZ2mod In on [0,27t]. Let f\ (t) be the puUback prob- 
ability density function of the nonuniform random variable Zi to [0,27t]. We let /^-(O = Prob(r < f) be the 
cumulative density function. Then 

Frit) = — / Prob(-/logZ2 e[s,s + f])/i {s)ds = — tfi {s)ds = t/{2%) (49) 

271 Jo 271 Jo 

Since Frit) ~ t /{2k), we see that T is uniform. Hence Z1Z2 is uniform. □ 
Proof of Theorem 1.7: First, let us checkthat {K(a) = 1} is t/a-measurable. To see this, note that the concurrence 
spectrum A-t(fl) may be expressed in terms of the coordinates tj as follows: 

dl = e^"«,di = e-'''-^''\dl = e2''2-2''i , • • -c/j = , . . . ,4_i = q-^^'n-i (50) 

Thus Lemma 111.2 induces a measurable condition on the tj. 

Continuing the proof, by direct calculation Z^ is a uniform random variable on the circle {zz = 1 } given that Z 
is such. Thus note that d^, d^, d\, • • • , are p — N /2 independent, uniform random variables by Lemma 111.4. 
It suffices to show that £+ \ = p independent, random variables on the circle have in their convex hull £^ 0°. 
Relabel d^ = Zq, c/j = Zi, • • • = Z^. 

Without loss of generality, say Zq ~ 1. Let C2 be the event that no Zi, Z2, • • • Ze is in the second quadrant 
{z=x + iy \ x<0,y> 0}, with C3 similar for the third quadrant {x < 0,y < 0}. Let D be the event that is in the 
convex hull of Zo, Zj, • • • , Zf . Then (NOT Ci n NOT C3) C D. Then Prob(NOT C2 n NOT C3) < Prob(D), and 

1 - Prob(D) < 1 - Prob(NOT Co and NOT C3) = Prob(C2 or C3) = (1 /2)^' (51) 
Hence as £ ^ 00, Prob(D) goes to 1. Hence the probability CH[X£.(\')] contains limits to 1. □ 



IV. CONCLUSIONS AND ONGOING WORK 

We have shown that there exists a generalized canonical decomposition of unitary operators on n qubits which 
may be used to study changes in the concurrence entanglement monotone. This decomposition closely resembles 
the older two-qubit decomposition when n is even, and it may be used to study the concurrence-entanglement 
capacity of generic unitary operators. The main result is that such a generic unitary operator is almost always 
perfectly entangling with respect to the concurrence monotone when the number of qubits is large and even. 

Ongoing work would attempt to extend the dynamical viewpoint taken in this paper. Specifically, the unitary 
operator describes the dynamics of a quantum data state, and the present techniques allow us to quantitatively study 
the dynamics of the concurrence entanglement measure. Similarly, we would wish to study the dynamics of this 
concurrence capacity of quantum computations in naturally defined families or sequences of such computations. 
As a separate topic, we might also study the failure of the concurrence function itself by quantifying how entangled 
a quantum state with zero concurrence may be. 

APPENDIX A: Computing the CCD When the Number of Qubits Is Even 

This appendix recalls how to compute the canonical decomposition in an even number n = 2p of qubits. Note 
that other arguments in the case « = 2 (12; 13) may be found in the literature, and that the present treatment is a 
straightforward genearlization of a matrix-oriented treament in the two-qubit case (14, App.A). It is included for 
completeness. 

The overall structure of the algorithm is contains two steps. 

1. Produce an algorithm for computing the decomposition SU (N) = SO{N) D SO{N) for D the diagonal sub- 
group of SU (N). We will refer to this decomposition as the unitary SVD decomposition henceforth. 

2. Recall £0 the standard entangler of Definition 11.9. Given av £ SU{N) for which we wish to compute the 
CCD, compute first the unitary SVD EqvEq = o\do2- Then we have a CCD given by 

V = {EQ0iEl){EQdEl){EQ02El) = hak2 (Al) 
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since ki = EqoiEq e K, k2 = Eq02Eq E K, and a = EqiJEq E A. 

Note that the unitary SVD decomposition exists due to KAK metadecomposition theorem, taking as inputs G = 
SU {N), 0Ai(^) = X, and a the diagonal subalgebra of 5u(A^). 

Before continuing to Step 1, we first prove a lemma. It is useful in computing particular instances of the unitary 

SVD. 

Lemma A.l For any p E SU {N) with p — p^, there is some o E SO{N) such that p — odo^ with d a diagonal, 
determinant one matrix. 

Proof: We first show the following. 

V a^b, symmetric real N xN matrices with ab = ba, there is some o E SO{N) such that oao^ and 
obo^ are diagonal. 

It suffices to construct a basis which is simultaneously a basis of eigenvectors for both a and b. Thus, say Vx is the 
A- eigenspace of b. For jc E Vx, b{ax) ~ a{bx) = Xax, i.e. xi-^ ax preserves the eigenspace. Now find eigenvectors 
for a restricted to V^' which remains symmetric. Thus we may find the desired o E SO{N), making choices of 
orderings and signs on an eigenbasis as appropriate for determinant one. 

Given the above, write p ^ a + ib. Now 1 = pp^ = pp = (a + ib){a — ib) ~ [a^ + b^) + i{ba — ab). Since 
the imaginary part of 1 is 0, we conclude that ab = ba. Hence a single o exists per the last paragraph which 
diagonalizes the real and imaginary parts. □ 

Suppose then that v = o\do2 is the unitary SVD of some v E SU {N). For convenience, we also label v = poj, the 
type AI Cartan decomposition ( •, thml.l.iii,pg.252) (2~, thm6.31.c). This is a generalized polar decomposition 
in which p — p^ , k E SO{N). Note that it is equivalent via Lemma A.l to compute v ~ pk, as the unitary SVD 
follows by V = {oidoJ)k = o\do2. Continuing to the algorithm for Step 1, 

• Compute p^ as follows: p^ = pp^ = po^o^p^ = vv^ . 

• Apply Lemma A. 1 to p^. Thus p^ ~o\d^o\ foroi E SO{N). 

• Choose square roots entry wise in d^ to form d. Be careful to ensure det d = I. 

• Compute p = oidoj . 

• Thus 03 = p^v, and v = p{oi) = oido^o^, = o\do2. 

This concludes the algorithm for computing the unitary SVD of Step 1. 

Step 2 is almost follows given the inline description. The reader may produce algorithms outputting Eq. 

Another question is computational efficiency. This is ongoing work, but we note immediately that an implemen- 
tation of the spectral theorem of Lemma A.l is required. This will be difficult with current technologies in 16+ 
qubits. Moreover, in the range of 50 to 60 qubits an even spread of the concurrence spectrum Xc (v) of Definition 
111.1 would make certain elements indistinguishable at 16-digit precision. 



APPENDIX B: Concurrence level sets and K orbits 

Mathematically, related measures are often easier to use than C„. For example, the concurrence quadratic form 
2n (l¥)) = C„(|v|/), |v|/)) with C„(|\|/)) = has smaller level sets than C„ itself. Moreover, it turns out that 

the normalized states within these level sets [2n] = llv) ^ Cn (Iv)) = ~} naturally orbits of the group 

K, which must then be false for C„. 

Suppose throughout n = 2p is an even number of qubits. For a vector v E C^, put Qai{v) = v^v, noting that 
Qc (Eov) ~ Qai{v). Moreover, for O E SO{N), we have the following: 

Qai (O ■ v) = Qc [EqOeI ■ (£0 v)] (B 1 ) 

Thus we may study level sets of Qai under SO{N) rather than study level sets of Qc under K. Now if v = vi + iv2 
is a decomposition into real and imaginary parts of a complex vector, note that Qx\{v\ + iv2) = v^v = (jvi |" — 
|v2|2)+2/(yi-V2). 
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Lemma B.l Label S^^ ' = {|\(/) ; (v(/|x)/) = 1}. We have the following orbit decompositions of the level sets 
Q^l (a) n S-^-^for any fixed a G C. 

7. Let t be real, and let v G Q^l (t) n S^^'K Then Q^l (f ) n S-^- ^ = [SO{N) ■ v]. 

2. Let a be complex, and let V e Qf^lia) nS'^'^-K Then Q]^l{a)r\S^^-^ = [SO{N)-v]. 

Proof: For the first item, write v = vi + iv2- Then vi • V2 is zero as a set of real vectors. Consider the subset of R^'^ 
given by |vi p — |v2p = f . Suppose now we have another pair of orthogonal vectors wi, W2 with |wi p — |w2p = f 
and 1^1^ + 1^21' = 1. Then \wi\^ = |vip = (1 -r)/2, thus |v2p = |w2p so that there is some O £ SO{N) with 

O Vi = Wi, (9 • V2 = W2- 

For the second item, suppose a = e'*? for some f G M. Now if v G Qai i^)' '^^en note that we have (2ai {e^"^^^v) = 
e^"^QAi{v) = e^"l'a = t. Conversely, if w G 2ai (0 we have e'^/^w G (a). Having established bijective phase 
maps between the two level sets, it must also be the case that the level set of a forms a single SO{N) orbit. □ 

Corollary B.2 The restricted action ofK to the normalized kets in any concurrence level set is transitive. Specifi- 
cally, suppose a G C, with |\|/) normalized with Qc ( |\|/) ) = OC. Then label S^'^^ ' = { ((|) |(|)) = 1 } the set of normalized 
kets. Per Equation Bl, we have K ■ |\|/) = Q'^^{a)r]S^^^\ 

We restate the result colloquailly. Should any two normalized states |(|)), |\|/) have the same concurrence, then 
there is some global phase e'^ so that |(|)) = e'^k\\\f) ioxkeK = Eq S0{N) E^. 

APPENDIX C: Concurrence is an entanglement monotone 

The «— tangle, defined to be T„(|x|/)) = C„(|x|/))^ has been proposed (1) as a measure of n qubit entanglement 
for n even. The «-tangle of a state \\\t), like the n-concurrence, assumes real values in the range < T„ < 1 and 
has been shown to be an entanglement monotone, meaning T„ is a convex function on states and is non-increasing 
under local operations and classical communication (LOCC). Most of our arguments focus on constructions more 
directly related to the concurrence C„ rather than the «-tangle T„ ~ [Cn)^- Therefore, for completeness, we show 
that the «-concurrence is, in fact, a good measure of entanglement. The monotonicity property of a function is 
established by considering its action on mixtures of quantum states encoded within Hermitian density matrices p 
with tr p = 1. See, e.g., (21). 

Definition C.l The «-concurrence can be defined on mixed states p using the convex roof extension: 



This minimization is over all pure state ensemble decompositions of the state p = L/t^/t|v'^) (V^^l- 

This definition is quite intricate. We point out the following remarkable result, not used in the sequel. 

Theorem C.2 (Uhlmann, (1; 25)) We may express C„(p) in closed form as follows: 



Here, the are the square roots of the eigenvalues (in non-increasing order) of the product pp where p = SpS ^ 

The necessary and sufficient conditions for a function on quantum states to be a entanglement monotone are 
delineated in (9). For the n-concurrence, they can be summarized as follows: 

• Qi > 0, and C„(p) = if p is fully separable. 

• C„ is a convex function, i.e. C„(ppi + (1 — p)p2) < pC„{pi) + {l — /?)C„(p2),V/7 G [0, 1] and pi,p2 Hermi- 
tian matrices of trace one 




(CI) 



C„(p) =max{0,A.o-A.i...-A.Af_i} 



(C2) 
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• C„ is non increasing under LOCC. Specifically, C„(p) > T^j PjCn{p j), where pj =AjpA'j/ pj are the states 
conditioned on the outcome j of a positive operator valued measurement (POVM) which occurs with proba- 
bility = fr[A]A^p]. 

Before proving this, we first establish the useful fact that the ^-concurrence is invariant under permutations of 
the qubits. Defining n„ to be the set of unitary operators corresponding to permutations on n qubits, we have; 

Proposition C.3 Forn even, C„(P|i(/)) = C„(|i|/)) VP G n„. 

Proof: Any permutation P on n elements can be written as a finite composition of transpositions on pairs of 
elements. Hence it suffices to show invariance under a single swap operation. Writing the swap operator between 
qubits j and k as 

1 / (g) + a-^ «) of + (g) + a^- (g) a? 
Sjk = (C3) 

we have for any state |\)/), 



= Msslm = c„(|v)). 



(C4) 



Here we have used the fact that Sji^ is real symmetric and unitary, and in the third equality we use the fact that 
[cr^ (g)cr^,a'ig)a[] =Ofora' G {cr*, a^,o'"}. This proposition necessarily implies that n+ C/T, where n+ is the set 
of unitary permutation matrices on n objects with +1 determinant, i.e. permutations composed of an even number 
of transpositions. □ 

Lemma C.4 C„(p) is an entanglement monotone. 

Sketch: For the first condition, one first checks that < C„ < 1 using the eigenvalue decomposition of the 
matrix S — (— /Oj)(— /o^,) • • • (— /oj,). Then any separable state can be realized by stochastic local unitaries acting 
on the fiducial separable state |0)„ = |0i . . .0„). Now, C„(|0)„) = and C„ is invariant under local unitaries per 
Proposition 11.1 and Theorem 1.3. To generalize from pure states to density matrices, recall Definition C.l. 
The second condition is shown by writing the minimal ensemble decompositions for p i and p2 separately as 

p min ^X,C„(|v')) + (l-/^) ^ min £p,C„(|^^)). (C5) 

These are not necessarily the minimal decompositions for the composite state p = pp\ + (1 — p)p2, therefore, 
C„(ppi + (1 -/:')P2) < pC„(pi) + (1 -p)C„(p2). 

Finally, we show that the «-concurrence is on average non-increasing under LOCC. First, because of permuta- 
tion symmetry of the concurrence we can consider operations on one particular qubit of the n qubit system, say 
the first. An arbitrary, trace perserving, completely positive map on a quantum system can written in the Krauss 
decompostion (26) as ^(p) = Jlj^jP^] where the positive Krauss operators satisfy the sum r ulei:^A]A; = l. The 
map can be composed of multiple operations with two operators at a time so we consider only two operators Aq 
and Ai acting on the first qubit. By the polar decomposition theorem, the operators can be written as Aj = Ujbj, 



where bj = yA^jAj is positive and m, is defined to be 1 on the kernel 3C of Aj and Aj|Aj| ' on !?C^. Physically, 

the map S{p) corresponds to a generalized measurement on p followed by a unitary operation conditioned on the 
measurement. Because of the sum rule, which corresponds to trace preservation, we can write Aq = Wo cos gX and 
Ai = wi singX for g G M and X a positive operator with unit trace. These operators are expressed in simpler form 
by diagonalizing X, viz. Aq = m„c/oV and Ai = uidiv, where Uj,\' G SU (2) and dj are real diagonal matrices with 

elements {q, r) and (^^1 — q^, vT~~r^). The average concurrence of a state p after the 2 outcome POVM is 

PoC„(po)+/5iC„(pi) = Pomin^;^^ |^t)j|^;^^j^i)(,|,*|^pIi,A.i:C„(Ao|v|/*^)/V/^) ^^^^ 
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States conditioned on the first outcome satisfy; 

C„(Ao|\|/)/V/^) ^ \{\\i\v'^d1^ulSuQdQv\\\i)\/ pQ 

= qrCn{v\^v))/po (C7) 
= qrC„{\\\t))/po, 

where in the second equality we have used Uq {—10^ )110 = —io^, and the third equality follows by invariance of 
concurrence under local unitaries. Similarly, C„ (A 1 1(|)) / = — q^){l — f'^)C„ ( |(|)) ) /p 1 • The result is, 

P0Cn{P0)+PlCniPl) = 'l^^'^^'^{l,,\^^|k)}\l^X,\^^|k){^^,k^^pLkhC„{\^\t'')) 

= (^r+V(l-?2)(l-r2))C„(p) 
< C„(p), 

with equality iff q^r, i.e. the Aj are stochastic unitaries. □ 
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